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In this paper, we propose a way to achieve protected universal computation in a neutral atom 
quantum computer subject to collective dephasing. Our proposal relies on the existence of a De- 
coherence Free Subspace (DFS), resulting from symmetry properties of the errors. After briefly 
describing the physical system and the error model considered, we show how to encode informa- 
tion into the DFS and build a complete set of safe universal gates. Finally, we provide numerical 
simulations for the fidelity of the different gates in the presence of time- dependent phase errors and 
discuss their performance and practical feasibility. 
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I. INTRODUCTION 



Within the last few years, quantum information has 
become one of the most promising and active fields in 
physics. In the commonly used model, a quantum com- 
puter consists of two-level systems, the qubits, in which 
information is stored in a binary fashion. Calculations on 
this information are achieved through the application of 
particular evolutions of the system, called quantum gates 
Thanks to quantum parallelism 0, which merely 
follows from the linearity of quantum mechanics, quan- 
tum computers are expected to be much more efficient 
than their classical analogues, in particular for simula- 
ting the behaviour of quantum systems Q and for solv- 
ing some "difficult" problems, such as factoring [4(. Un- 
fortunately, none of the various experimental proposals, 
comprising NMR Cavity Quantum Electrodynamics 
@, and trapped ion Q implementations, has succeeded 
in fulfilling all the requirements one has to check in order 
to design a valuable quantum computer, known as Di- 
Vincenzo's criteria ||. In particular, decoherence which 
arises from the interaction of the system with its envi- 
ronment, remains a major obstacle to the feasibility of 
quantum computing. 

Quantum gates may be implemented with neutral 
atoms using either short range collisions or long range 
dipolc-dipolc interactions. The proposal of Jaksch et al. 
_9| suggested using the strong dipole-dipole interactions 
of highly excited Rydberg atoms for fast quantum gates. 
In a recent article [T(j, one of the authors provided a 
detailed study of the Rydberg scheme using optically 
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trapped 87 Rb atoms. In particular, it was shown that 
the errors due to the trap setup itself can be made quite 
small by a proper choice of physical parameters, allowing 
for fast and reliable single- and two-qubit gates. 

In this paper, we elaborate on this proposal by con- 
sidering the effect of collective random dephasing errors, 
which for instance stem from the uncontrolled action of 
exterior fields. In Sec. UH we briefly present the physi- 
cal implementation proposed in [lfj, as well as the er- 
ror model we choose to address and we show how to 
protect information through encoding into a Decoher- 
ence Free Subspace (DFS), the existence of which merely 
follows from the symmetry properties of the errors. In 
Sec. IIII1 we build a complete set of universal protected 
gates, which allows us to perform any computation with- 
out leaving the DFS. In Sec. HVi we provide numerical 
simulations for the fidelity of the different gates in the 
presence of time-dependent phase errors. The practical 
feasibility of these gates is then discussed; in particular, 
the limitation arising from spontaneous emission is ad- 
dressed. Finally, in Sec. [Vj we give our conclusions and 
the perspectives of our work. 



II. A NEUTRAL ATOM QUANTUM 
COMPUTER 



The quantum computing proposal we shall consider 
throughout this paper has been put forward recently in 
[To| . The physical qubit consists of a 87 Rb atom re- 
stricted to the hyperfine states 

|0) = |5S 1/2 ,F=l,m F = 0) 
|1) = |5S 1/2 ,F = 2,m F = 0) 

which form the logical basis (see Fig. []}. After pre- 
cooling in a magneto-optical trap (MOT), the different 
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III. PROTECTION AGAINST THE ERRORS 
THROUGH ENCODING INTO A DFS 



F=2 



F= 1 



ID 



A. Secure storage of the information in a DFS 



-2-1012 



10) 



-1 o 1 



FIG. 1: The qubit basis states. 



qubits which constitute the computer are captured in an 
egg box-style potential created by a Far-Off-Resonance 
Trap (FORT) (for a detailed presentation of the techni- 
cal aspects of the physical apparatus, see 

In this setting, the single qubit rotations are imple- 
mented through Raman- like transitions between |0) and 
|1) via an off-resonance excited state, whereas the two- 
qubit phase gate is performed according to Jaksch et al.'s 
proposal Q which relies on the large dipole-dipole inter- 
action between Rydberg atoms. Combining single-qubit 
rotations and two-qubit phase gates leads to the standard 
set of universal quantum gates, comprising the single- 
qubit Hadamard (H) and tt/8 (T) gates, as well as the 
two-qubit CNOT gate : 
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Analysing the error sources due to the trap setup 
(background gas collisions, scattering of the trapping 
light, heating due to laser noise, ...), the authors of [1G | 
show that, given a proper choice of the physical param- 
eters, single- and two-qubit gates can be performed at 
MHz rates with decoherence probability and fidelity er- 
rors at the level of 10 -3 . 

In this paper, we shall consider the effect of random 
phase errors which affect different qubits according to the 
same (a priori time-dependent) error Hamiltonian 



E(t) = h 



eo(t) 

£l (t) 



Typically, this kind of error model stands for the action 
of parasitic external fields which induce uncontrolled and 
unwanted energy shifts. If no protection scheme is used, 
the unknown differential phase shift induced by E (t) on 
the qubit states rapidly leads to a complete loss of coher- 



To protect information from quantum errors, one can 
resort to active schemes such as quantum codes[l|, the 
most famous example being the stabilizer codes |11| . It is 
also possible to take advantage of the symmetry proper- 
ties of the interaction between the system and its environ- 
ment in order to passively protect information from the 
effects of decoherence; this is the basic principle of the 
Decoherence Free Subspace {DFS) strategy [12||, which 
consists in encoding information into subspaces immune 
to errors. The explicit construction of DFS has alr eady 
been achieved for certain collective error processes [13| . 
and even experimentally implemented in quantum optics 
[lij and trapped ion [15[ setups. Moreover, universal 
computation within these DFSs has been shown to be 
possible from a theoretical point of view 16]. Here, we 
shall deal with the practical implementation of such a 
DFS in the neutral atom system described in [Tol ] . 

In the case of collective dephasing, a DFS can be 
straightforwardly identified in the Hilbert space of a two- 
physical-qubit system. Indeed, as the two qubit states 
{\0dfs) = 1 01) , \1dfs) = 1 10)}- are affected in the same 
way by phase errors 



DFS ) = |01) e 
t 



1 



DFS) 



|10) 



/ '(eo(r)+ei(r))dr 
e - i /o( e o(T)+ei(T))<ir 1 10) 



the subspace they span is clearly left invariant by the 
errors. This shows that a qubit of information can be 
safely stored in the state of a two-physical-qubit system. 

From a practical point of view, in order to protect the 
information contained in a physical qubit initially in the 
state \ip) — Co |0) + c\ |1), one simply adds a second auxil- 
iary qubit, initially prepared in the state |1) and performs 
a CNOT gate, which achieves the DFS encoding 



|V) |1) = c |01)+ Cl |11) CN S T c |01)+ Cl |10) = |V 



'DFS/ 



Decoding is straightforwardly performed through apply- 
ing the CNOT gate again. Moreover, a protected N- 
logical-qubit memory can be constructed just by associat- 
ing N protected cells, which requires 2N physical qubits. 

To process the information safely stored in the DFS 
we now need to design a new set of universal gates, com- 
prising the single- logical-qubit Hadamard and ir/8 gates, 
as well as the two-logical-qubit CNOT gate. This point 
is actually not obvious. Indeed, the gate implementa- 
tions proposed in [ToT | cannot be used any longer as they 
would make the system leave the DFS, resulting in erro- 
neous calculations. We thus have to resort to new phys- 
ical primitives. In the following we review two secure 
processes which constantly remain in the DFS. We then 
show how to combine them to implement universal pro- 
tected computation. 
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FIG. 2: Level scheme for the dephasing gate. 



B. Physical primitives for universal protected 
computation 

The first primitive consists of a phase gate, obtained 
through applying a laser to the first atom of the logical 
qubit. The Rabi frequency is denoted and the laser 
frequency LOd is detuned by the quantity from the 
frequency ui of the transition |1) — > |e) where |e) denotes 
an excited state (A^ = tod — ui, see Fig. [5]). After a 
time td — fp 5 -) wnere &>R,d = vNdF + Al an d n is a 
positive integer, the state 1 10} picks up the phase factor 
e ltp , where ip = nix (l + n^)j while the state 1 01 > is 
left unchanged. In the DFS basis, this transformation is 
represented by the matrix 



P(<p) = 



1 

e 1 ' 



a pair of atoms. The frequencies Ui are assumed slightly 
detuned from the transitions \i) — > |r), i — 0,1, respec- 
tively, where \r) is a Rydberg state of the atom. The 
detunings are denoted (see Fig. Oi)) 

A = luq — <^>0r and A' = ujq — u>\ — cJoi 

where u>i r is the frequency of the transition \i) — > |r). Let 
us emphasize that the two atoms involved in the trans- 
formation are not bound to belong to the same logical 
qubit, i.e. they do not have to be in a superposition of 
the states |01) and |10). 

As is well known [17| . Rydberg atoms exhibit huge 
dipole moments which lead to large dipole-dipole inter- 
actions responsible for dipole blockade @. We shall as- 
sume that the dipole-dipole interaction is only signifi- 
cant when both atoms are in their Rydberg states |r). 
In other words, we shall assume Vdd = hA rr \rr) (rr\. 
When the fields are applied to a pair of atoms, transi- 
tions may occur between the different resonant physical 
qubit product states, in agreement with energy conser- 
vation. We observe that the states 1 00) and |11) only 
couple to themselves to second order, and, as the dou- 
bly excited Rydberg state \rr) is shifted far off resonance 
due to the dipole-dipole interaction, the only two paths 
actually coupling 1 01} and 1 10} are the two fourth-order 
'M' systems drawn on Fig. [3p). 

Assuming |fi„| , \^\ < |A| , |A'| , |A' - A| , \A rr \, one 
can extract the effective dynamics of the subspace 
1 1 00) , 1 01) , 1 10) , |11)} through a perturbative approach 
[l8| which yields the evolution operator, expressed in the 
interaction picture, 

U eff (t, n 0: fii, A, A', A rr ) = e -*( A o-A')* 



Given an arbitrary phase ipo one can choose the different 
physical parameters such that (p — ipo and \Ad\ ^> Nd|, 
which ensures that the excited state |e) remains essen- 
tially unpopulated during the process. 

The second primitive is implemented by applying two 

laser fields Ei (t) = \ ^Eq^ 1 ^-^ + ccj , i = 1, 2, to 
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FIG. 3: a) Level scheme and laser couplings for one atom, b) 
The two 'M' paths coupling |01> and 1 10) . 



It is important to note that if the two atoms involved are 
initially prepared in a DFS state, i.e. in a superposition 
of 1 01) and |10), the transformation U e ff leaves them in 
the DFS and the associated gate 



R{6) 



cos (9) i sin (9) 
i sin (9) cos (9) 



with 9 — 



2 



is thus safe. But it is also interesting to note that if 
the two atoms considered are not in a DFS state, as for 
example in the state |00), they pick up a phase factor, 
which can be controlled via the Rabi frequencies and de- 
tunings. This observation will be used in the following 
to implement a two-logical-qubit phase gate. 

To conclude, let us emphasize that the transforma- 
tions we have just considered are not completely error- 
free. They are safe in the sense that they do not make 
the system leave the DFS, the information thus being 
constantly protected from the effect of the errors. But 
the errors will of course affect the primitives P and U e f / 
(and thus the gates we shall build by combining them) 
through the detunings A, A', A rr which will take slightly 
different values from the ideally expected ones. However, 
the detunings being large, the effect of reasonably small 
errors on the system will be quite unnoticeable. The 
limitation of the fidelity for the quantum gates will be 
numerically analyzed in the following section. 



C. Construction of the universal safe quantum 
gates 

Let us now see how the previous primitives can be 
used to implement universal safe computation. The 
single-logical-qubit ir/8 (T) and Hadamard (H) gates are 
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FIG. 4: The U e ff gate acting on atoms 1 and 3. 



straightforwardly obtained according to 
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We now consider two logical qubits, consisting of the 
physical qubits (1,2) and (3,4) respectively. The two- 
logical-qubit DFS basis is 

|00) DFS ^ |0101) 
\01) DFS = |0110) 
\W) DFS = |1001) 



111) 



DFS 



|1010) 



One might think that in order to implement a logical 
two-qubit gate a four-atom process must be involved. It 
turns out, however, that a two-atom operation is suf- 
ficient to implement a controlled phase gate between 
two logical qubits. Thus consider the transformation 
(t, fio, Hi) A, A', A rr ) involving atoms (1) and (3) 



(1,3) 
eff 



from different DFS pairs, cf. Fig. [4] This operation 
leaves the subspace {|00) DFS , |01) DFS , |10) DFS , |11) DF5 } 
invariant and is given by 
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in the DFS basis. If we adjust the different parameters 
in such a way that 
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FIG. 5: The CNOT DFS gate. 



corresponding to the conditional phase gate in the DFS 
subspace. It is now straightforward to implement the 
protected CNOT gate 

CNOTdfs ={I® H dfs ) • U phas ^ DFS ■ (I (8 #dfs) 

cf. Fig. [5l which completes our set of protected universal 
quantum gates. 



IV. NUMERICAL SIMULATIONS 

As already pointed out, the gates we have just imple- 
mented are not error-free, as their fidelities are limited by 
the errors which modify the detunings involved. Sponta- 
neous emission from the excited states |e) and |r) also re- 
stricts the performances of our gates, and has to be taken 
into account. In this section, we numerically estimate the 
influence on the fidelity of random phase errors and spon- 
taneous emission and compare the protected gates to a 
set of unprotected gates. Finally we discuss the practical 
interest and feasibility of our proposal. 



A. Performances of the protected gates 

In this subsection, we first describe how phase errors 
and spontaneous emission are modelled in our system. 
Then we introduce the fidelity measure we shall use 
throughout this section. Finally we explain how simu- 
lations are performed, and discuss our results in detail. 

Here we assume phase errors on |1), |e) and |r) and 
model their effect by the single-atom Hamiltonian 

E(t) = h[e 1 (t)\l){l\+e e (t)\e){e\+e r (t)\r){r\] 

Moreover, we suppose that phase errors can be described 
by an Ornstein-Uhlenbeck process for which 19]: 

e(t + dt) = e(t) - -e(t)dt + y/cG(t)Vdt 

T 

where r is the relaxation time, c is the diffusion constant 
and G(t) is the unit Gaussian variable. The Ornstein- 
Uhlenbeck process is characterized by the correlation 
function 

(e(t)e(t')) = ^e-\ t - t '\/r 



in steady state. We assume that the errors for the various 
levels are correlated, which seems reasonable for back- 
ground magnetic field perturbations. Thus e\(t) is gener- 
ated using the approach described above while e e (t) and 
e r (t) are found from e e (t) = a e t\{t) and e r (t) — a r e\{t). 

Spontaneous emission from the excited state (used for 
the implementation of the P gate) and the Rydberg 
state can be taken into account by adding the term 
— ih {^f\e){e\ + ^|r)(r|) to the single-atom Hamiltonian. 
For the numerical implementation we have considered 
\e) = \r). The radiative line width of the Rydberg state, 
2ir, is estimated to be smaller than 10 kHz for n > 65 

An appropriate fidelity measure for the performance 
of the gates in the presence of stochastic errors is calcu- 
lated in the following way. For an N x N density matrix, 

p, the time evolution is given by p — > U pU^ . Repre- 
senting the matrix p as a vector p this can be rewritten 

such that ~p — > M~p with M being an TV 2 x N 2 matrix. 
To each instance of t\ (t) (and correspondingly e e (t) and 
e r (t)) corresponds a certain M. In order to obtain an 
average over the various phase error configurations, we 
can now simply average over the corresponding Ms and 
subsequently calculate the fidelity from: 

where p ldea i = U ldeal pllJ deal , M av denotes the average 
over M and p = In practice, we evaluate F for 

a large number of pure states \ip), and in the figures we 
present the minimum value of F over these input states. 

As shown in the previous section, the different univer- 
sal gates can be obtained by combination of the two prim- 
itives P (tp) and U e ff (t, fioj A, A', A rr ): it is thus 
necessary to perform these primitives with high precision. 
The expression obtained for P (if) is exact: it implies 
that, in absence of errors and spontaneous emission, if we 

tune {n, Qd, A^} such that if = rnr ( 1 H — , ^ d I is 

\ vl^l +|a<*I / 
the desired phase factor, we shall exactly get the expected 
phase gate P '(</>)■ On the contrary, the expression we 
obtained for U e ff is only perturbative. It means that if 
we choose the physical parameters {t, Qq, f2i, A, A', A rr } 
so that the effective parameters Eq.JT]|3]) take the de- 
sired values (to implement a conditional phase gate or 
an R(8) gate), the exponential of the full-Hamiltonian 
will be slightly different from the desired gate, even in 
the ideal case. Before dealing with errors and sponta- 
neous emission, we refine the parameters by a numerical 
search in the neighborhood of our first guess provided by 
the analytical perturbative approach. Once we have a 
faithful gate in the ideal case, we keep the same set of 
parameters and include the new terms discussed above 
in the Hamiltonian in order to run our simulations. 

Fig. [6] and [7] show the performances of the different 
gates when spontaneous emission is not/is included, re- 
spectively. The fidelity is plotted as a function of 
which corresponds to the variance of the detuning due 
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FIG. 6: Performance of the universal set of gates without 
spontaneous emission. Dash-dotted: Tdfs, dashed: Hdfs, 
solid line: U p hase,DFS , dotted: C NOTdfs ■ Parameters for 
the gates: P gate: Q. d /2n = 5 MHz, A d /2-K pa -49.81 
MHz, t(7r/4) pa 1.00 fj,s. 7?(f) gate: Q /2tv pa 3.91 MHz, 
fh/27r w 1.97 MHz, A/2tt pa 60.34 MHz, A'/27r pa 30.70 
MHz, A rr /2n = 100 MHz, £ pa 120.20 fis. U phase , DF s gate: 
n /27r pa 3.93 MHz, fii/27r pa 1.96 MHz, A/2tt ps 60.48 MHz, 
A'/27r pa 30.05 MHz, A rr /2iv = 100 MHz, t pa 938.62 (js. For 
the generation of the phase errors the parameters are: r = 1 

jUS, Qe = Or = 1.5. 
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FIG. 7: Performance of the gates including spontaneous emis- 
sion with 7r/27r = 5 kHz. Parameters and notations are the 
same as for Fig. [6] 



to perturbations. As is evident in Fig. [7| spontaneous 
emission severel y li mits the fidelity. 

Note that in [lOj the transitions to the Rydberg states 
are implemented by going off- resonantly via \5P\/ 2 ) ■ This 
worsens the problem of spontaneous emission, as the 
R and U p hase,DFS gates are thus quite slow and the 
radiative linewidth for the \5P\n) state is quite large 
(j 5 p 1/2 /2tt = 5.7 MHz). It is thus necessary to be 
strongly detuned from the level \bP\/ 2 ), but as this re- 
duces the coupling to the Rydberg levels, the intensity of 
the Raman beams should at the same time be increased. 

A further source of decoherence for the experimental 
setup proposed in [lfj is the motion of atoms in the traps: 
as the interatomic distance, i?, varies so does the dipole- 



dipole interaction, which scales as 1/R 3 . Based on esti- 
mates provided in [l(| a variation of 20 % for the value of 
A rr is expected. We performed a numerical simulation 
assuming that A rr is harmonically varying which demon- 
strated a dramatic reduction of the fidelity for realistic 
parameters. Atomic motion in the traps thus poses a 
serious limitation for the gates. To improve the perfor- 
mances of the gates the atoms might be cooled further or 
the distance between the traps could be increased, which 
would, however, reduce A rr . 



B. Comparison with unprotected gates 

In order to assess the performance of the universal set 
of DFS gates suggested in this paper, we compare them 
to a universal set of unprotected gates. 

If we remove the restriction that the system should 
remain in the DFS at all times, a controlled-phase gate 
between two atoms can be carried out as suggested in 
Q in the regime of a large interaction strength. Com- 
bining this gate with two single-atom Hadamard gates, 
which are obtained through single-atom rotations be- 
tween states |0) and |1) via an excited state, results in a 
two-atom CNOT gate. A universal set of gates is then 
given by: 

H DFS = CNOT 1>2 ■ (H (g> I) ■ CNOT 1>2 
C NOT DFS = CNOT h3 ■ CNOT 1A 

where CNOTi.j is a CNOT gate with atom i being the 
control qubit and atom j being the target qubit. The 
Tdfs gate is simply implemented in the same way as for 
the protected gates. 

In Fig. [5] we plot the performance of the unprotected 
Hadamard and CNOT gates. As expected, these gates 
are much more sensitive to phase errors than the pro- 
tected gates. Spontaneous emission from the Rydberg 
states also severely limits the fidelity for the unprotected 
gates, as the controlled-phase gate involves resonant tran- 
sitions to the Rydberg levels. 

The durations of the gates are, however, smaller than 
for the protected gates. Thus, with the parameters for 
Fig. [8] the durations of the gates are t(HDFs) = 15.5 fj,s 
and t{CNOT DF s) = 14 MS- 

Another advantage of the unprotected gates as op- 
posed to the protected ones we have suggested is that 
they are less sensitive to variations in A rr , and thus to 
atomic motion. 



V. CONCLUSION 

In this paper we have identified a set of logical qubit 
basis states for a neutral atom decoherence free subspace 
and a corresponding set of protected universal gates. Nu- 
merical simulations demonstrate that the proposed set of 
gates is much more robust against phase errors than a 
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FIG. 8: Performance for the unprotected gates. Dashed: 
H DFSl solid line: CNOT DFS . Parameters: Q/2tv = 0.5 MHz, 
A rr /2n = 100 MHz, r = 1 fjs, a e = a r = 1.5. 



for computations, as information is then protected dur- 
ing storage of the qubits. 

The work presented in this paper naturally leads to 
further developments. First, we would like to study re- 
source optimised DFS encoding. Here, in order to build 
a reliable N — 2-qubit register we have merely associated 
two two-physical-qubit protected cells, but it is clear that 
for an iV-qubit quantum memory [N > 3) the number of 
physical qubits required by this straightforward scheme, 
i.e. 2N physical qubits, is much larger than actually 
needed. This new construction of a DFS will induce new 
problems, as regards practical encoding and processing 
of the information stored in the DFS, the resolution of 
which will probably require the use of quantum control 
techniques, such as nonholonomic control (20l |. Moreover, 
we also plan to extend our results to more general error 
models as for instance position-dependent errors. 



set of unprotected gates. They are, however, also much 
more affected by the motion of the atoms in the traps. 
Therefore, one has to consider what the worst source of 
decoherence is for a given physical situation and subse- 
quently choose which set of gates to use. Even if the 
set of protected gates is assessed to be unfavorable, one 
should still consider using the logical qubit bases states 
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